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Abstract 

The lifetime of localized surface plasmon plays an important role in many aspects of plasmonics 
and its applications. In small metal nanostructures, the dominant mechanism restricting plasmon 
lifetime is size-dependent Landau damping. We performed quantum-mechanical calculations of 
Landau damping of a surface plasmon in a metal nanoshell. In contrast to the conventional picture 
based on the electron surface scattering, we found that the damping rate vanishes in the limit of 
thin nanoshell. The physical origin of this behavior is the reduced magnitude of bright surface 
plasmon electric field inside the metal shell. We also found that, due to electron scattering on two 
metal surfaces, the damping rate exhibits pronounced quantum heats with changing shell thickness. 
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I. INTRODUCTION 



Lifetime of localized surface plasmon (SP) in metal nanostructures is one of fundamental 
problems in plasmonics that has been continuously addressed for more than 40 years 
The importance of this issue stems from one of the major objectives of plasmonics - genera- 
tion of extremely strong local fields at the nanoscale. The range of physical phenomena and 
applications related to this goal cuts across physics, chemistry, biology, and device applica- 
tions. A small sample of examples includes SP-enhanced spectroscopies of molecules near 
metal nanostructures, such as surface-enhanced Raman scattering (SERS) and fluorescence 
resonance energy transfer (FRET) and their numerous biomedical applications ^ , focusing 
of extremely large electric fields in metal nanoparticle (NP) clusters ^, and quantum coher- 
ent SP generation in nanostructures (SPASER) j^. High Ohmic losses in bulk metal due to 
strong electron-phonon interactions impose limitation on the quantum yield of metal-based 
plasmonic devices, which can, to some extend, be remedied by reducing the metal compo- 
nent size. However, at the lengthscale below ~ 10 nm, new limitations on SP lifetime and, 
consequently, on quantum yield arise due to quantum-size effects. Among those, the most 
important is the Landau damping (LD) of SP - decay of SP into the Fermi sea electron-hole 
pair. This process is made possible by relaxation of the momentum conservation due to 
electron surface scattering. The electron-hole momentum uncertainty, 6p ~ h/a, where a is 
the characteristic size of nanostructure, translates into the corresponding energy relaxation 
scale, 6E ~ p6p/m, m being electron mass in metal, and this leads to the following estimate 
of SP lifetime, 

r~A^, (1) 

a 

where vp is the electron Fermi velocity (hereafter we set h = 1). The first explicit quantum- 
mechanical calculation of P for spherical NP by Kawabata and Kubo [1] confirmed the 
estimate ([1]) with a replaced by particle radius R and, for NP in vacuum, yielded the value 
A = 3/4. While the subsequent studies indicated that the coefficient A varies in the range 
(0.5-2.0) depending on dielectric constant and chemical composition of surrounding medium, 
the size dependence ([T]) was thoroughly confirmed experimentally for spherical particles Q]- 
The corresponding quality factor, 

Q = ^, (2) 
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where uisp is SP frequency, decreases linearly with R. This causes a rapid reduction of, e.g., 
SERS enhancement factor, M oc Q^, while the other quantum-size effects, such as electron 
density spillover or surface screening, play a subdominant role j^. 

The above qualitative arguments have been subsequently used for estimates of SP lifetime 
in nanostructures of various shapes, notably in metal nanoshells The single- electron 
picture of SP damping, based solely on surface scattering, inevitably leads to a SP damping 
of the form ([1]), with a replaced by the smallest spatial scale of a confined nanostructure. For 
example, in a nanoshell with shell thickness d, the characteristic electron scattering length 
is ~ d, and Eq. ([T]) predicts a very short SP lifetime for thin nanoshells. First attempts to 
measure SP lifetime in nanoshells seemed to be consistent with Eq. ([T]) with estimated values 



A = 2.0 



and A = 0.5 . However, more recent systematic studies of optical spectra 



of single gold nanoshells indicated no size-dependence at all {A = 0) for the extracted SP 



lifetime 



io|. 



There is also a physical argument that renders Eq. ([T]) invalid for nanostructures of general 
shape. Indeed, LD is a many-body phenomenon, rather than a single-particle one, and it 
involves the SP own electric field by which an electron-pair is excited. Therefore, the SP 
lifetime must be highly sensitive to the field magnitude inside the NP. Note that for a solid 
sphere, the SP electric field magnitude is a frequency-dependent constant that is independent 
of NP radius. This is the reason why Eq. ([1]) holds remarkably well for spherical NP in a very 
wide range of their radii, from several tens of nm down to nanometer-sized small clusters 
In contrast, for other NP shapes, the magnitude of SP electric field inside the NP can 
change with NP size. As we show below, the latter effect can, in fact, be much stronger 
than surface scattering contribution, resulting in a drastically different, from Eq. ([1]), SP 
damping rate. 

In this paper, we calculate the LD of bright SP modes in a metal nanoshell with dielectric 
core. We demostrate that in the limit of thin nanoshell, the LD rate vanishes linearly with 
the nanoshell thickness d. Furthermore, for small overall nanoshell sizes, SP lifetime exhibits 
quantum beats as a function of shell thickness, caused by the interference between electron 
scattering amplitudes near inner and outer nanoshell boundaries. We also show that despite 
the redshift of the SP frequency, the nanoshell quality factor increases with decreasing d. 
The latter finding indicates that, with appropriate choice of geometry, the quantum-size 
limitations imposed on plasmonics can be sidelined, and the quantum yield of even small 
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plasmonic devices can be made limited only by the Ohmic losses. 

In Section [IT] we obtain a formal expression for LD rate in terms of plasmon eigenmodes in 
metal nanoshell with dielectric core. In Section UTTl we evaluate of nanoshell internal energy, 
and in Section IIVI we evaluate the dissipated power by calculating the imaginary part of 
electron polarization operator in semiclassical approximation. The analysis and numerical 
results are presented in Section |Vl and Section |Vl] concludes the paper. 



II. PLASMON EIGENMODES IN A NANOSHELL 



We consider metal nanoshell with inner and outer radii Ri and R2, respectively, in a 
medium with dielectric constant ea- In long-wave approximation, the local potential $(r) is 



determined from self-consistent equation [llj 

<l>(r) = V5(r) + j dr,dr2U{r - r^)P{r,, r2)$(r2), (3) 

where (p oc r^Y^Mi^) is an external potential {Ylm{^) are spherical harmonics), P = P'+iP" 
is polarization operator, and U{r) is Coulomb potential. We are interested in plasmonic 
eigenmodes that are described by Eq. ([3]) without external potential, which we write as 

(l + 47rx)$ = 0, (4) 

where x is a susceptibility operator defined as 

X$(r) = e^ j dridr2G{r - ri)P(ri, r2)$(r2), (5) 

and G{r) = — f/(r)/47re^ = —l/Anr. For optical frequencies, the real part of polarization 
operator can be approximated by its local limit and decomposed into the core, metal and 
outside dielectric contributions, P' = P'^ + + P^, with 

i^'(r,r') = e-2V[x.(r)V5(r-r')], (6) 

where Xc{r) = XcO{Ri - r), Xm{r) = XmO{r - Ri)9{R2 - r), Xd{r) = XdO{r - R2) are local 
susceptibilities, Xi = i^i ~ l)/47r. 

In the absence of damping, plasmon energies are obtained by setting P" = and solving 
Eq. (jl]) in the basis of eigenfunctions of composite operator x' = GP': x'Flm = ^lFlm- 
The eigenvalues are found by setting 

1 + A-kXl =em + B± VB^ + C = 0, (7) 



where 

R — -^^cm - {L + l)emd ^ _ L{L + 1) 2L+1X ^ /q\ 

^" 2(2L + 1) ' ^"(2LTiF^ 

K = R1/R2 is the shell aspect ratio, and we denote Sais = Sa — ^p- The eigenf unctions 
Flm = -^L('")^LM(r) are obtained from continuity of -Fi(r) across all interfaces 

Note that Fq satisfy orthogonality condition (/i|P'|z/) = S^uP^; /i is composite index (LMj), 
where j = it in Eq. (JTj) is plasmon type (bright or dark, respectively). 

Plasmon decay into electron-hole pairs is described by incorporating the imaginary part 
of polarization operator into Eq. x = GP = GP' + iGP", as follows. First, Eq. (jl]) is 
multiplied from left by {fi\P' and, using that P'x = x'Py we obtain 

(1 + 47rA^) (/i|P'|$) + AmX^{^i\P"\^) = 0. (10) 

Then, expanding the potential over eigenstates of x', ^ = J2 ^l^lmi S^^ 

(1 + 4vrAi)itCi + 4vrai i^^'C^ = 0, (11) 

L'j' 

where 

P'l^^ = ^ E (^A^JI^'I^'M'J') (12) 

MM' 

is the reduced, due to spherical symmetry, matrix element of P". Admixture of plasmon 
modes by P" being weak, the non-diagonal terms in Eq. flTTl) can be neglected. The new 
resonance condition takes the form 1 + 47r(A-^ + (5A-^) = 0, where 5A]^ = iX-'j^P'l" / P'l is an 
imaginary correction to the eigenvalue. Its appearance implies a shift of eigenmode frequency 
into the complex plane, Xl{!^l + if^i/'^) = Xl{ijJl) + i{VL/2) [dXiiuJi) / duj^]. Equating the 
last term to 5Xl-, we obtain the following expression for the decay rate 

_ 2XlP'U^l) _ Ql 
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where we suppressed bright /dark indices. Here Ql is the dissipated power of plasmon 
eigenmode, 

Ql = -ujlP'll = J^A^ ■ Elm, (14) 

Elm = —^Flm and ^lm being the plasmon electric field and the current, respectively, and 

is the plasmon eigenmode energy stored in a metal-dielectric composite. 

III. CALCULATION OF EIGENMODE ENERGY 

Consider first the palsmon eigenmode energy, . The electric field integral in Eq. (fT5|) 
reduces to boundary contribution, 

j dr [e{r) - 1] Ei^, = -47r [/??FL(i?i)4'^ + RlFLm^f^] , (16) 

where aj*'* are surface charges at ith interface. Those are given by cr^^^ = ecm£^L(-Ri)/47rec 
and cr^^'* = emd^L(-R2)/4vrerf, where SiiRi) = —dFi{Ri) / dRi are the interfacial electric fields 
(on metal side). The prefactor in Eq. (fT5l) can be split as 



-ALULide'JduL), (17) 



UJL_d\L_ 

Xl duL 
where 

Al = -AZ' (9AL/9e„) = 47r (SAi/Se^) . (18) 

For metal shell, we use the Drude dielectric function Sm = Sdb — ^p/^iyJ + ^7), where 
Edb is the inderband dielectric function due to c?-band to sp-band transitions, Up is bulk 
plasmon frequency, and 7 is the bulk damping rate. Then we have ojL^de'j^/dujL) = 
2{uji/ujp)^ ^iuji/ujj^"^ + (7/co'p)^] ^ ^ liupl uji)'^ for optical frequencies {ujl ^ 7). Putting 
all together, we finally obtain 



Ul = — ^ ( — 



2 



L + 1 



Sl{R2)Rt + -^Sl{Ri)R{ 



(19) 



where we used £l{Ri) = -L (eJemRi) Fl{Ri) and £l{R2) = {L + l) {ed/£mR2) Fl{R2) 
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The dependence of dipole SP energy on aspect ratio is plotted in Fig. [T] for gold nanoshell 
in vacuum, and for gold nanoshell with AU2S core in water. Importantly, with increasing 
core radius, the stored energy remains comparable to (vacuum) or exceeds (dielectric) that in 
a solid particlle. This can be understood by noting that a nanoshell is analogous to a (high- 
frequency) capacitor formed by two concentric spherical surfaces. For bright dipole mode, 
the opposite charges are located in different hemispheres, so the energy is proportional to 
the core-shell particle volume; a dielectric core then leads to an increase of the stored energy. 
In contrast, for dark modes, the opposite charges are placed at inner and outer boundaries, 
so the energy vanishes as with the volume between capacitor plates decreases. 

Simple expressions can be obtained for a metal shell in vacuum, = = 1- The 
corresponding eigenvalues are given by 

Al = AlIb^ - l)/47r, Al = ^[it VTT4L{L + T)^ / {2L + 1)] , (20) 



yielding the well known resonance frequencies ujl/^p = y [1 + AL{edh — 1)]~^- In this 
case, the integration in Eq. ( fTSll is effectively restricted to metal region and we recover the 
known expression for time-averaged internal energy in dispersive media 



Ul = f^p^ [ drEl,, = ^ldr 
Svr oul J on 



duJis'jujL, r) 
dujr. 



Km (21) 



where we used V ■ {sElm) = 0. Then Eq. (1191) can be evaluated using eigenf unctions (E 
with r]L = Lk^ [L + 1 - (2L + 1)^^]"^ The resuh reads 

' ^ L{L + l)/€2^+i 



1 + 



[L + 1 - (2L + l)ALf 



(22) 



For thin shells, d/i?2 = (1 - k) < 1, we have {ul/uJpY ^ Al ^ {I - k)L{L + 1)/(2L + 1) 
(for bright plasmon mode), so that Ul ~ R2{'^L + l)^/47r(L + 1)^, while for solid particle 
we have Uf^ = R2{Ledb + L + l)/47r, yielding Ui/Uf ~ 0.6 [note that Ul vanishes in the d.c. 
limit Lj/j <1 corresponding to d/R2 < {•y/UpY]. 

IV. CALCULATION OF DISSIPATED POWER 

We now turn to the calculation of dissipated power, Ql = —^lP'll- Within RPA for 
imaginary part of electron polarization operator, it is given by 

Ql = ^^Y. mFLM\a')\'5{e^-e^,+u:L). (23) 

Maa' 



7 



where the sums run over occupied (a) and unoccupied {a') electronic states with wave- 
functions ^/'a(r) = ipni{r)Yim{r) and energies e„i {n is radial quantum number), and the 
factor 2 accounts for the spin degeneracy. The angular part factorizes out and Eq. (1251) can 
be written as 

Ql = 2nujL J2 <^w\M!:iyi,?5{tni - tn'l + ujl), (24) 

nn'lV 

where M^;^,^, = {nl\FL\n'l') and afi, = ^ | J dfiFiM^m^i'm' |^ are radial and angular 

Mmm' 

contribution, respectively. The latter is non-zero only for I = I' ±L, and for typical /, /' ^ L, 
can be approximated as a^, ~ 5;/'//27r. 

Let us first qualitatively discuss the LD of the L = 1 SP in a solid particle of radius R. 
In a finite system, the dipole matrix element, {nl\r\n'l') = ^^(n/|p|n7') is non-zero because 
the momentum is not a good quantum number. For low-energy transitions across the Fermi 
level {ul <^ Ep) we have {nl\p\n'l') ~ pp/{n — n') ~ eopp/ui, where eo is separation between 
levels at the Fermi level with the same (or fixed difference of) / (eo = hvp/R for a solid 
particle). We then obtain {nl\p\n'l')'^ ~ [e^ppY /irP'ujf. The sums over n, ra' in Eq. flM|) are 
transformed into energy integral over a strip of width uji around the Fermi level, resulting 
in a factor ujipf, where pi ~ Cg ^ is the partial density of states (DOS) at the Fermi level. 
Finally, the sums over angular momenta is dominated by the maximal / in a quantum well 
of width i?, l^ax ~ VfR) contributing another factor of {ppR^. Putting all together, we 
obtain Qi ~ R"^ [Ep/ujif' . Note that the level spacing eo cancels out, and dissipated power is 
determined by the nanoparticle surface area. At the same time, the energy of dipole plasmon 
(with Fi = r) is proportional to the particle volume, Wi ~ (i?^/e^) (tUp/tUi)^, yielding 
F ~ [Ep/upf (e^/-R) ^ Vp/R. We emphasize that the 1/R dependence of F originates from 
the distribution of surface charges rather than from quantum-size corrections. The dissipated 
power is simply proportional to the surface area available for electron scattering, but it is 
insensitive to the structure of electron energy spectrum, as indicated by the cancellation of 
electronic DOS in the transition matrix elements. A similar cancellation takes place also 
for a shell geometry so that Ql is also proportional to its surface area. On the other hand, 
the bright plasmon mode energy Ul is proportional to the core-shell particle volume, so 
one would expect a similar size-dependence of F as for a solid particle of the same overall 
size. However, as we show below, LD damping in thin nanoshells is even weaker to due the 
enhancement of screening, caused by the SP frequency redshift, that pushes the bright SP 
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electric field out of the metal shell region. 

To evaluate Ql, we represent nanoshell confining potential as 3D quantum well with hard 
boundaries at -Ri and R2, V{r) = Vq6{Ri — r)9{r — R2). The matrix element of plasmon 
potential F^m is dominated by surface contribution, {a\FLMW) = -^{ci \H [H^Flm\\ \a') ~ 
^^^(alVF^A/ ■ Vl^|a'), and, after separating out angular part, takes the form M^i^^,i, = 
^ [ipniiRi)ipnn'{Ri)£LiRi) - i'niiR2)ipnn'{R2)£L{R2)]- For dccp cuough wcU, matching of 
wave-functions across the well boundaries gives y/2mVoilJni{Ri) ~ ~'^nii^i) (here prime 
stands for derivative), and the matrix element takes the form 

M^^,, = [i^URlWn'Vm^Lm - i^UR2Wn'Vm^L{R2)] . (25) 

L 

The two terms describe excitation, by the SP electric field, of a Fermi sea electron-hole pair 
with the momentum transfer to inner and outer boundaries. Correspondingly, Ql can be 
decomposed as Ql = Q]^ + ~ 2Q^^, where 

■£L(i?.)^L(i?,) 5Z V';z(^*X'i(^^X;(^,)^:.';(^,)^(en/ - e^'i + Ul), (26) 



^ Inn' 



and we used that a^, ^ Sii/l/27!-. For typical electron energies e„i ~ Ep, we can use 

R2 

semiclassical approximation for electron wave-functions, ipniir) = ^/Am/piTi sin J pidr, 

r 

where pi{e,r) = A/2me — (/ + l/2)^/r^, and Ti{e) is the period of classical motion be- 
tween two turning points. Consider first the outer surface contribution, Qf^. In this case 
'ijj'^i{R2) = — ^/ Ampi{R2) / Ti. When plasmon energy is larger than fixed-/ level spacing in 
the nanoshell, eo = vp/d, the sums in Eq. (!26l) can be replaced by integrals, J2 ^ I d^Pii.^) 

n 

(with e < Ep,(' > Ep), where the partial DOS pi{e) = dn/deni is related to the classical 
period as pi = ti/ 271. The result reads 

Ep 

Qf = 4^y2^ [ depi{e,R2)pi{e + u;L,R2). (27) 

Note that the partial DOS cancels out and Q^^ depends only on electron and hole velocities 
at the outer surface R2. In the energy integral, the integration variable is first shifted as 
e Ep + e — ujl/2, where e now changes in the interval (—^7^/2,6^^/2), and then rescaled 
to X = e/ujp. The sum over I is replaced by integral restricted by maximal value I ~ PfR2 
that is determined by the condition p;(e, i?2) > 0. After rescaling to s = l"^ /{ppR2)'^, it 
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contributes a factor proportional to the outer surface area. The result reads 

Qf = ^Sl{R,)g{u^/Ep), (28) 



1/2 



where g{^) = 2 J dx J dsf{^, x, s), with f{^, x, s) = y (1 + — s) — ^^/4, is dimension- 

-1/2 

less function of order 1 normalized to g{0) = 1. 

Turning to the inner surface term, Q^, the main contribution into the r.h.s. of 
Eq. (1261) comes from the terms where the turning point lies outside of potential well, 
i.e., pi{e,Ri) > (otherwise ipni{Ri) are exponentially small). In this case, we have 



il)'^i{Ri) = —{—l)'^\jAmpi{Ri)/Ti, where the sign factor (—1)" accounts for the parity of 
electron wave-function with n — 1 nodes between Ri and R2. The rest of the calculation is 
carried in a similar way, and the result, 

QY = ^SliRMuj/Ep), (29) 

is proportional to the inner surface area. 

Consider now the interference term Q]^. Using the above boundary values of ip'^i{Ri), we 
write 

A£LiR,)£L{R,) ^ /(-!)»-"' 

Ql = > , zJ^l[(^nh en'l>0{enl - Cn'i + ^^l), (30) 

^ riieniMen'i) 

where J^z(e„i, e„//) = ^/piieni, Ri)pi{en'i, Ri)pi{enu R2)pi{^n'h ^2)- The alternating sign re- 
flects the parity of excited electron-hole state. As ujl changes (e.g., with changing aspect 
ratio), the relative parity of electron and hole states separated by energy ujl changes too, 
leading to a different sequence of alternating signs in the sum in Eq. ( l30l) . This leads to 
oscillatory behavior of Q}^ . The number of states contributing into the sum on Eq. fl30|) 

being large, the oscillations are relatively weak and can be accounted, in the continuum 

r 1 

limit, by a replacement (— 1)"~" = cos7r(n — n') = cos vr J depi{e) . Then takes the 
form 

Ep r e+uii^ 

QF = ^#44^E^ / de:Fde,e + u,)cos vr / de'pi{e') , (31) 

Ef~uil L e 

where / is restricted by the condition pi{e,Ri) > 0. After shifting integration variables as 
e — >• Ep + e — ujl/2 and e' Ep + e + e' and rescaling, Eq. fl5Ul) takes the form 

Ql = ^^Sl{Ri)Sl{R2)G{ujl/Ep), (32) 
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where 



1/2 1/2 



-1/2 



-1/2 



, (33) 



R2 

and the partial DOS is given by p^e) = f / -f^r = (l/Svre) [i^sP^e, ^2) - i?iPKe, Ri)] with 
I = y/sippRi)- For ijjjEp <^ 1, the x'-integrals is easily evaluated, yielding 

1/2 i-f^x 

G{i) = 1 y rfx y c^sa/ /(^, X, X, K^s) COS X, s)u;l/ eo] , (34) 

-1/2 



where w;(^,x,s) = a/ /(^, x, k^s) -ka/ /(^, x, s) with/(^,x, s) ^ 1+^x-s (here eo = t;F/^2) 

1/2 

Rescaling s by 1 + ^x, Eq. (IMI) factorizes as ^(0 = / (ix (1 + ^x)^''^ S'(^, x), where 

-1/2 



5'(^,x) = 2 / ds\/ (1 — s)(l — K^s) cos a(^, x) ( Vl — /t^s — /tVl — s 



(35) 



with shorthand notation a(^,x) = iiLii^ltQ)^! + ^x. With substitution s = 1 ^sinh a, 



5" is brought to the form 



4(1 - K 



2\2 



Q!0 



da (sinh a cosh a) cos a(^, x) Vl — /t^e " 



(36) 



where sinhao = k/\/1 — k^. For a(^,x) ^ 1, the integral is dominated by the upper limit, 
and for thin shells, 1 — k ^ 1, corresponding to ao > 1, can be evaluated as 

sin [a(l — k)] 



sin(aVl - K^e"""^ 



-4- 



a(l — 



(37) 



With the above 5* and after change of variable t = v^l + ^x, the expression for G{^) takes 
the form 

t+ 

G(0 = -^ / (38) 



where t± = a/1 ±^/2, and we introduced D = {1 — K)ujL/eo = ULd/vp. Note that even 
though for ^ <^ 1 the integration interval is small, the integrand is still an oscillating function 
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since D ^ 1 and so the product can be arbitrary. In this straightforward 
evaluation yields 

G(0 ^ -45l^5l|^. (39) 

Putting all together, we finally obtain 

Ql = [Sim + ^'Slm - 2^'Sl{Ri)Sl{R2)G] . (40) 

The last term in the above expression oscillates as a function of n. Oscillations are caused 
by the change of electron-hole state parity as its energy changes with d, and reveal them- 
selves through the interference between electron scatterings from outer and inner surfaces. 
The oscillations period depends weakly on shell thickness throught ujl{k,), and their am- 
plidude slowly dies out with increasing d. Note that although the interference correction, 
Eq. (!30|) . was obtained for thin shells, the full dissipated power Q^, Eq. (HOj) remains a good 
approximation in the entire range of k, since for small k the interference term is also small. 

V. LANDAU DAMPING IN NANOSHELLS 

We now turn to SP LD rate, = Ql/Kl, with Ql and Ul given by Eqs. 1^ and (|40D. 
To simplify the analysis, consider the case of a metal shell in vacuum [ec = Sd = !)• Then 
we obtain 

1 + K'^ql - 2n''qLG 



-e2 



(41) 



where qi = £l{Ri) / £l{R2) is the ratio of electric fields at the boundaries. 



Lk^ -{L + l)riL K^-^A 



K[L-{L + l)7]Ln^+^] Az.-(1-k2^+i)(L + 1)/(2L + 1)- ^^^^ 
The expression f l4T|) is our central result. The SP LD rate is proportional to the energy of 
Fermi sea electron- hole pair created by SP (second factor), multiplied by the phase space 
of the transition (first factor). The last factor describes the relative contribution of the 
nanoshell inner surface and includes the interference correction. Importantly, the electron- 
hole Coulomb potential is screened by metal dielectric function at SP frequency. 

For solid nanoparticle (k = 0), the plasmon eigenfrequency is determined from Lem{ujL) + 
L + 1 = 0, and we recover the result of Kawabata and Kubo, 

rl -mV^- (43) 



TTLU^ J R2 4 i?2 
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Note that although the r.h.s. resembles the level spacing (at constant /) in spatially-confined 
system, this resemblance is coincidental because Eq. fH3|) contains no information about 
electron energy spectrum. 

For thin nanoshell, the situation is drastically different. The redshift of plasmon energy, 
ul implies large (negative) value of 61(001), which leads to overscreening of the electron-hole 
Coulomb energy. As a result, F vanishes in the thin shell limit. The precise behavior of 
Tl can be deduced from (1 — k) <^ 1 asymptotics, qi ~ —L/{L + 1) and ~ ~ 
-(1 - k)-\2L + 1)/L{L + 1), yielding 

i.e., F^, decreases linearly with shell thickness. 

In Figs. [2] and [3] we show the calculated LD for dipole SP as well as the corresponding 
quality factors in the entire range of aspect ratios for Au nanoshells in vacuum with outer 
radii 30 nm and 10 nm, respectively. The decrease of F is accompanied by the increase of 
Q despite the SP frequency redshift for thin nanoshells. The asymptotic behavior of quality 
factor, Q oc c?"^/^, follows from the weaker dependence of the nanoshell SP frequency, 
oc \/d. For small d, both F and Q exhibit oscillations as a function of d, with period 



UJ 



sp 

di = vpTi, where Tl is the SP period. The latter itself weakly depends on d, resulting 
in a slow increase of oscillation period with d, however undetectable due to the reduction 
of their amplitude. The oscillations are quite pronounced for smaller overall nanoshell size 
and could be observable for typical experimental range of aspect ratios (0.6-0.8). Note that 
these oscillations should be distinguished from those observed in solid NP [l^ ; in the latter 
case, the oscillations originate from size-quantization of electron energy levels in a confined 
nanostructure, while here they are due to quantum interference of electron wave-functions 
at different nanoshell boundaries. 



VI. CONLUSION 



In summary, we performed the first quantum-mechanical calculation of surface plasmon 
Landau damping in a metal nanoshell, and found that it vanishes in the thin shell limit. 
The physical origin of a long bright SP lifetime is that its electric field magnitude inside a 
thin metal shell is reduced. Obviously, such a situation is not limited to nanoshells with 
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spherical symmetry considered here, but can take place in other structures too. Therefore, 
the quantum-size limitations on plasmon lifetime in small metal nanostructures are, in fact, 
spurious and can be avoided by an appropriate choice of geometry. 
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FIG. 1: SP energy in Au nanoshell in vacuum (solid line) and in Au2S/Au core-shell particle in 
water (dashed line) relative to solid nanoparticle of the same overall size. 
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R,/R, 

FIG. 2: Calculated SP LD rate (a) and quality factor (b) for Au nanoshell with outer radius 
R2 = 30.0 nm in vacuum relative to solid Au NP of the same overall size. Solid horizontal line of 
unit height is added for comparison with solid NP. 
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FIG. 3: Same as in Fig. [21 but with R2 = 10.0 nm. 
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